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Abstract We investigate a new symmetry of the large deviation function of certain time-integrated 
currents in non-equilibrium systems. The symmetry is similar to the well-known Gallavotti-Cohen- 
Evans-Morriss-symmctry for the entropy production, but it concerns a different functional of the 
stochatic trajectory. The symmetry can be found in a restricted class of Markov jump processes, where 
the network of microscopic transitions has a particular structure and the transition rates satisfy cer- 
tain constraints. We provide three physical examples, where time-integrated observables display such 
a symmetry. Moreover, we argue that the origin of the symmetry can be traced back to time-reversal 
if stochastic trajectories are grouped appropriately. 
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1 Introduction 



In Nature most processes are not in thermodynamic equilibrium. For example, whenever a system 
is exposed to a flux of matter or energy in the stationary state, then it is generally not possible to 
describe it with standard methods of equilibrium statistical mechanics. For systems out of equilibrium 
the probability distribution of configurations is normally not known and a theory which allows one to 
calculate macroscopic quantities is not available. Therefore, general results valid for nonequilibrium 
processes are of great theoretical interest. This applies in particular to fluctuation relations [TUTU], 
which are fairly general statements valid for any system out of equilibrium. 

A large variety of nonequilibrium systems can be modelled as continuous-time Markov jump pro- 
cesses, meaning that the system jumps spontaneously from one classical configuration to the other 
with certain rates. If these rates obey detailed balance, the system relaxes into a equilibrium state. 
However, if detailed balance is broken, the system will relax into a nonequilibrium stationary state 
where the probability current between microstates is non-zero [TT]. Maintaining these non- vanishing 
probability currents requires an external drive which continuously produces entropy in the environ- 
ment. Remarkably, this entropy production can be quantified without knowing the explicit structure 
of the environment [12] . It turns out that the average entropy production is zero if and only if detailed 
balance is fulfilled, meaning that entropy production can be used as an indicator of nonequilibrium. 

Since the entropy production depends on the specific sequence of microscopic transitions (the 
so-called stochastic path), it is a fluctuating quantity with a certain probability distribution |13j . 
A fluctuation relation is an equation that restricts the functional form of this distribution. There 
are two different kinds of fluctuation relations, namely, finite time fluctuation relations and infinite 
time fluctuation relations (see [5]). Finite time fluctuation relations describe the distribution of the 
total entropy (system + environment) and hold exactly for any time interval [13J . Relations of this 
kind include the Jarzynski equality |14j and Crooks relation |15j . On the other hand, infinite time 
fluctuation relations are asymptotically valid for the entropy produced in the environment (which is 
also known as action functional) [5J. Here we will deal with this infinite time relation, and we refer to 
it as the fluctuation theorem or the Gallavotti-Cohen-Evans-Morriss (GCEM) symmetry. 

Large deviation theory |16H19| is the appropriate mathematical framework to investigate the fluc- 
tuation theorem. Using these methods the GCEM symmetry can be recast as a symmetry of the large 
deviation function for the probability distribution of the entropy. This symmetry does not yet allow 
us to calculate macroscopic observables, but at present it is the most general result for systems out of 
equilibrium. 

Even though the fluctuation theorem is very general, one might argue that it is also very specific 
in the sense that it is valid only for one particular functional of the stochastic path, namely, the 
entropy produced in the environment. It is therefore interesting to find out if there are other physically 
relevant functionals with a similar symmetry. As a first step in this direction, we recently demonstrated 
that the height of an interface in a certain growth model defines a physically relevant time-integrated 
current different from the entropy with a symmetric large deviation function [20] . Interestingly, this 
symmetry can only be observed for a particular system size of the model because only then the network 
of microscopic transitions acquires a particular form. 

The objective of the present paper is to generalize the new symmetry found in |20j . We prove that 
for a class of jump processes with a particular network of states we can find time-integrated currents 
different from entropy displaying a symmetric large deviation function. This symmetry is similar to the 
GCEM symmetry because it restricts the form of a large deviation function related to a time- integrated 
current. However, it is important to note that our symmetry is different from the GCEM symmetry 
because it refers to a different time-integrated current and has a slightly different physical origin. 

It is known that the GCEM symmetry is a direct consequence of the fact that the entropy is given 
by the weight of a stochastic path divided by the weight of the time-reversed path. In this sense, the 
origin of the GCEM symmetry is related to time-reversal. An interesting question, which we address 
here for a particular case, would be to explain the origin of the symmetry. We show that it is also 
associated with time-reversal, but in a more hidden way: it comes to light only when we perform an 
appropriate grouping of stochastic trajectories. 
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The organization of the paper is as follows. In Sec. 2 we define time-integrated currents and briefly 
review the fluctuation theorem. In Sec. 3 we prove the new symmetry relation. Some physical examples 
where our symmetry appears in physically meaningful time-integrated currents are presented in Sec. 4. 
Before concluding we discuss the origin of the symmetry for a simple four states system in Sec. 5. 



2 Time-integrated currents and the fluctuation theorem 



Continuous-time Markov jump processes are defined by a space of microscopic configurations c e fl in 
which the system evolves by spontaneous transitions c — > d at rate w c ^ c i. The probability P(c, t) of 
finding such a system at time t in the configuration c evolves according to the master equation 



p( C , t) = j2 ( p ( c '. *R'^ C - p{c, t) Wc ^ c ,y (i) 



d_ 

dt 

c'^c 



For simplicity we assume that stationary probability distribution exists P(c) = lim t _ i . 00 P(c,t). If this 
distribution obeys the condition of detailed balance P(c)w c ^ c ' = P(c')w c i-y c the stationary state is an 
equilibrium state, otherwise the system is out of equilibrium. 

A stochastic trajectory during the time interval [to,tf\ is a sequence of M jumps 

^M,t ■ c(t ) -> cfa) -> c(t a ) c(t M ) (2) 

taking place at times t\,t2,- ■ • , ijif G [to,t/]. Note that, the length of the time interval T = tf — to 
is given while the number of jumps M is a random variable that may assume different values for 
different trajectories. In what follows we assume that all microscopic transitions are such that, i.e. 
w c ^f C i =/= <^> Wc'^-c 7^ 0, meaning that any stochastic path can be reversed. 

Time-integrated currents 

A time-integrated current is a functional of the stochastic trajectory 

Af-l 

J0M,t] = ^2 0c(U)^c(t i+1 ) , (3) 
i=0 

which changes its value by 9 C -> C > whenever a jump from c — > c' occurs. The increments are assumed to 
be antisymmetric, i.e. C _> C ' = —6 c i^ c . Using the master equation the expectation value is 



(J)= / ' dtJ2^c'P(c, 



t)w c ^ c > . (4) 



In the stationary state this expression reduces to (J) = I 1 ^]^, ^ c /P(c)w MC », i.e. the current in- 
creases on average linearly with T. Since the system relaxes towards a stationary state, in the limit of 
T — > oo the quotient J/T tend to a constant, which is 

lim — y^6 c -^. c ,P(c)w c -> c >, (5) 

c.c' 

More specifically, one expects that the corresponding probability distribution P (if = x) becomes more 
and more peaked around this value as T —> oo. Assuming that the large deviation principle holds, the 
large deviation function of this probability distribution is defined by 



r limP^=^ =exp[-TI(x)}. (6) 
The function I(x) measures the rate at which the current deviates from its average value. 
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Entropy production and fluctuation theorem 

A very prominent time-integrated current is the entropy. If we have a jump process describing a physical 
system in contact with external reservoirs, this quantity describes the amount of entropy which is 
generated by the external driving in a fictitious external environment [H1EI1HI]- More specifically, 
each transition c — > d changes the entropy in the environment by In ™ C J* C ' . Therefore, the accumulated 

entropy production along a stochastic path ~C M,t is a time-integrated current of the form 

M-l 

J s 0M, t }= V In jgsgihiSgaO . (7) 

The fluctuation theorem reads [5] 

I a (x) -Is(-x) = -x, (8) 

where I s (x) is the large deviation function associated with limj^oo P{^f = x )- This property of the 
entropic current is also referred to as the GCEM symmetry. 



Determining the large deviation function 

The master equation (fTJ) can be rewritten in the form 

j t P(c,t) = -J2£cc>P(c',t), (9) 

c' 

where C is the Markov generator with elements 

(io) 

Here 

A(c) = (11) 

denotes the escape rate from configuration c. For each time- integrated current of the form ([3]) , one can 
now define a modified generator C(z) by 

_/ --uv^ c exp(-z0 c /^ c ) ifc^c' , , 

L[z)cd - \ A(c) if c = c' • 

The scaled cumulant generating function I(z), of the respective time- integrated current J, is defined 

by 

lim (exp(-zJ)) = cxp(-T/(z)). (13) 

T->oo 

It can be shown that I(z) is given by the minimum eigenvalue of the modified generator (fT2f [5]. 
Moreover, the Gratner-Ellis theorem [T6HT8] states that I(x) is given by the Legendre-Fenchel transform 
of I(z), i.e., 

I(x) = max z (^I(z) - xz^j , (14) 

with z real. Note that 1(0) = because in this case C(z) reduces to the Markov generator C with the 
minimum eigenvalue 0. The GCEM symmetry (jSJ in terms of the scaled cumulant generating function 
of the entropy I s (z) is 

I s (z) = I s (l-z). (15) 

The advantage of dealing with the scaled cumulant generating function is that it is easier to calculate 
in several situations. In the present case it corresponds to determine the minimum eigenvalue of the 
Perron- Frobenius matrix (fT12j) . Note that, the curve I s (z) has a convex shape, vanishes at z = and 
z = 1, and reaches its minimum at z = 1/2. 
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Cm c n Q-i 

Fig. 1 Artificial network of configurations c\ with periodic boundary conditions. The process can jump between 
configurations that are connected by a line. 

Non-entropic time-integrated currents 

Is it possible to find other time-integrated currents with a GCEM-like symmetry? Previous works (see 
e.g. [9j[22]) have shown that such currents do exist. However, these examples are unsatisfactory in so 
far as the proposed currents differ from the entropy only initially while they become proportional to 
the entropy in the long time limit. If we call such a current J r , this means that I r (z) = I r (E — z). 
Moreover, since the current J r becomes proportional to entropy in the long time limit, I r (zE) = I s {z), 
i.e., the rescaled scaled cumulant generating functions of J r and J s have the same functional form. 

The main objective of the present paper is to show that it is possible to find symmetric currents 
which differ from the entropy even in the limit T — > oo so that their rescaled scaled cumulant generating 
functions differ from I s (z), even tough they are both symmetric and touch the horizontal axis at the 
same points (see below in Figs. [4] [6l and [8]). 



Counting the degrees of freedom 

Before proceeding let us point out that there are certain restrictions that reduce the degrees of freedom 
in the space of time-integrated currents. As such currents are specified by an antisymmetric matrix 
6 cc i — —B C ' C there are in principle N(N — l)/2 degrees of freedom, where N is the number of states. 
However, not all of them arc independent. To see this let us consider the elementary timc-intcgratcd 
currents 

A/-1 

Jc-yc'0M,t] = ^2 (8c,c(ti)8c',c(U +1 ) ~ $c' ,c(U)$c,c(t i+1 )) (16) 

from which all other currents can be constructed by linear combination. The current J c -> c i [C M,t] is 
simply the number of transitions c — > c' minus the number of reverse transitions d — > c along the 
stochastic path T) M,t- The sum over all destinations 

N c = ]T J c -,c> (17) 

c' 

is just the number how often the system reaches the configuration c minus the number how often this 
configuration is left, hence N c can only take the values and ±1. This implies that 

d 

in the limit T — ¥ oo, meaning that these particular linear combinations of elementary currents do 
not contribute to the large-deviation function. This reduces the degrees of freedom by the number of 
independent relations (|18p . which is maximally N. 
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3 A symmetric time-integrated current different from entropy 



In this section wc prove that for jump processes with a particular network of states and suitably chosen 
transition rates one can define a current with a symmetric large deviation function which differs from 
the one for the entropy. The structure of this network is shown in Fig. [1] It consist of configurations 
c l k organized in columns labeled by a lower index k = 0, . . . , Q — 1, each of them including rife different 
configurations labeled by an upper index i = l...n&. Spontaneous jumps are allowed only between 
configurations in neighboring columns with periodic boundary conditions, as indicated by straight 
lines in the figure. Moreover, we assume that the number of columns Q is even and that even columns 
carry only a single configuration, i.e. no = n-2 = n± = ... = 1 . This forces the system to go through 
periodically arranged bottlenecks of single configurations. 

On this network of configurations we consider the current 



k=0 



n k iifc + i 

E E 

i k =i i k+1 =i 



j 



(19) 



where 9k are numbers and J c _>. c ' are the elementary currents defined in Eq. (|16p . This current increases 
by 9k if the process process jumps from column k to k + 1 and decreases by the same amount if the 
system jumps in the opposite direction, no matter which of the configurations within a column is 
selected. 



Counting cycles 

Assume that the process starts at a particular configuration, say cq. Depending on the yet unspecified 
transition rates, the process will perform a random walk from column to column. Whenever it returns 
to its starting point cq, it is easy to see that the current defined above will take the value J r = mO, 
where O = Yl^o ® k an< ^ m £ 2. The number m tells us how often the system completed a cycle 
through all columns 0—^1— > 2 — y ... — > Q — 1 — >• 0. Therefore, if the rates are chosen in such a way 
that the random walk through the columns is biased to the right, m will be on average positive. With 
this picture in mind it is intuitively clear that the expectation value of m, in the long-time limit is 
related to the average current through each of the bottlenecks. 



To prove this intuitive argument, we apply the restriction (|18[) to each of the configurations in the 
network. For single but multiply connected configurations at even columns this restriction tells us that 
the incoming current is equal to the outgoing current in the long time limit, i.e. for even k we have 

E 15 ^ = E^r ±i - (20) 

For the simply connected configurations in columns with odd k, we have instead 

(21) 



T T 

where i — 1, . . . , n^. From the above relations it is easy to show that the current in the large deviation 
regime is given by 

Y = eE%^. (22) 

i=l 

This result proves that the large deviation properties of the current do not depend individually on 
the contributions 9 k but only on their sum = ^2 k ® k - This means that all currents of the form (|19D 
are proportional in the long-time limit and therefore characterized (up to rescaling) by the same large 
deviation function. 



G 





Fig. 2 Local transition cycle (left) and complete transition cycle through all columns (right). 



Structure of the characteristic polynomial 

We now prove that for suitably chosen transition rates the scaled cumulant generating function of J r 
exhibits the symmetry (|T5|) . Following Ref. [5] we first consider the characteristic polynomial 



P(z,x) = det(xl - t(zj) = ^ sgn(?r) ]~[ (xS c . /rTc - Cjz) c ^ c ) 



(23) 



where the sum runs over all permutations tt of the available configurations. Here C(z) denotes the 
modified generator (|12[) with the matrix elements 



t(z) d . c 



exp(-z 6 c ^ d ) + 5 c ^\(c) 



(24) 



where A(c) is again the escape rate Since the off-diagonal entries represent the possible microscopic 
transitions only those permutations will contribute to the determinant which correspond to a set of non- 
intersecting transition cycles. In the network of configurations shown in Fig. [T] there are three types of 
closed transition cycles, namely local cycles which do not change the current J r , and cycles extending 
over the whole system in positive or negative direction, changing the current by ±<9 (see Fig. [2]). 
In the determinant this means that permutations corresponding to local cycles are z-independent 
since the exponential factors drop out. Conversely, complete cycles extending over the whole system 
contribute to the sum with terms which are proportional to e ±z0 , respectively. We can therefore split 
the characteristic polynomial (|2"3"]l into three parts 



P(z,x) = f(x)e-* e + f(x)e*» + g(x) 



(25) 



Labelling each complete cycle cq — > cy 1 —> c 2 
the first two functions can be expressed as 



Cq_1 cq by a multiindex i := (ii, *3, . - . , *q-i) 



fix) = Y, T ^ 



ax) = Y, T i R &) 



(26) 



where Ti and T\ are the products of all rates along the cycle in forward and backward direction, 
respectively. Explicitly, 



T{ = w 



Ti = w 



'Q-i w *Q-i. 



(27) 

Moreover, R comprises all diagonal entries coming from the configurations that are not involved in the 
cycle: 



Q/2 n 2k -i 

w*) - n n 

fc=l 32k-i=l,jafc-i^»2fc-i 



-A(4-i 1 ) + ^) 



(28) 



There is also a complicated expression for g(x) but, as we will see below, this function is not needed 
for finding the symmetry. 
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Symmetry condition and derivation of the constraints 



Clearly, a symmetry of the characteristic polynomial P(z, x) = P{E — z, x) is a sufficient condition for 
the symmetry I r (z) = I r (E — z) of the minimal eigenvalue. Because of Eq. ([25]) this means that the 
condition 

/(*) = exp(E9)J(x) (29) 

implies a GCEM symmetry of the current defined in (|19p . Since f(x) and f(x) are polynomials in x of 
order N — Q, where N = E& n k is the total number of configurations, we can compare the coefficients 
on both sides. Equating the leading order x N ~® we obtain 

cxp(Ee) = |^p- . (30) 

l^i 1 i 

The above formula says that exp(EO) is given by the sum of the product of the transition rates of 
each possible forward cycle that goes trough all the columns (increasing the current by 0) divided by 
the sum of the product of the transition rates of each possible backward cycle. 

The comparison of the N — Q remaining orders yields a set of constraints for the transition rates. 
To this end we rewrite Eq. (|28p as a power series 

N—Q 

Mx)= ^(-l) n ^- Q - n a„(A0, (31) 

n=0 

where a n (Yi,Y 2 , Y N -q) = Y.i<i 1< i 2< ... < i n <N~Q Y h Y h ■■■ Y l n is tne elementary symmetric polyno- 
mial and where 

M = { x (4 k ) fc = l,3...,Q-l; j k = l,...,n k ; Jk £ i k } (32) 

is a set of N — Q arguments, consisting of all escape rates which are not part of the cycle labeled by i. 
Inserting this power series into Eq. (|2"r?)) and comparing the coefficients in Eq. ([2^]) one is led to N — Q 
constraints for the transition rates of the form 

EiTi = EiTitJniA) n=l2 ,...,N-Q (33) 



Simple solutions of the constraints 

Even though these constraints appear to be complicated, they arc fulfilled trivially by setting 

A(4) := X k (34) 

for all k = 0, . . . , Q — 1 and i = 1, . . . , n k , meaning that all configurations in the same column have the 
same escape rate. In this case the function Ri(x) = R(x) does no longer depend on the specific choice 
of the cycle labeled by the multiindcx i, meaning that the polynomials f(x) and f(x) reduce to 

f(x)=R(x)J2 T ii J(x) = R(x)J2Ti. (35) 

i i 

Clearly these functions satisfy equation (|29[) , which implies in a symmetric characteristic polynomial 
and therewith a GCEM- like symmetry of the current (|T9)) . 

Another trivial solution that fulfills the constraints is the following. Let us take a column k that 
has more than one state and define the quantity 

W i k W i k 

Ft" = Cfc ^ fc Cfc k+1 , (36) 

^ W ^ k W ^ k 
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where ik = 1, 2, . . . , n^. If this quantity is constant in all columns with more than one site, i.e. 



F k =Fl k for i k = l,...,n k , (37) 

then the ratio 

rp Q/ 2 

=- = n^ fc -i (38) 

1 i fe=l 

is independent on the cycle labeled by the multiindex i because the reversal any complete cycle changes 
the product of the rates always by the same factor, hence we arrive at the symmetry condition 

/Q/2 \ _ 

/(*)= I; f[ i^fe-i J f(x). (39) 

Finally, we get a larger class of solutions by mixing the conditions (|3~4"|) and (|3T|) . i.e., some of the 
Q columns with more than one state have constant escape rates while the remaining columns have 
constant Fk- Using similar arguments one can again show that the ratio /// is constant and, therefore, 
the constraints ((33)) arc fulfilled. Nevertheless, we note that not all the solutions of the constraints 
equations are of this type. 



The relation between J r and J s 

Depending on the transition rates the current J r may be proportional to entropy. If this is the case, 
then J r displays the GCEM symmetry. However, when J r is not proportional to entropy in the large 
deviation regime and still with a symmetric large deviation function, then we have a symmetry different 
from the GCEM symmetry In the following we derive the condition on the transition rates such that 
J s is proportional to J r . 

For the network of states shown in Fig. [T]the entropy current is given by 

fc=l,3,5,... \i=l c fc -»c fc _i j=1 c k+1 ^c k J 

Using relation ([2~Tj). in the long time limit the above term divided by T becomes 

'^= r E E^^ln Wl - (41) 

fe=l,3,S,... i=l c fc + 1 ->c fc c fc ^c fc _! 

Let us first consider the second solution of the constraints (|3T|) . In this case the quotient of rates 
appearing in the logarithm does not depend on i, therefore, using relation (|20[) we obtain 

£ = ^( E ^fE^V (42) 



,fe=l,3.5, 



In this case J r is proportional to the entropy and thus exhibits a GCEM symmetry. However, we can 
also conclude in opposite direction that for systems obeying the constraint equations (|33|) in such a 
way that (|37l) is not satisfied for all k, then the current J r is not proportional to J s . In this case we 
expect to have a new type of symmetry different from the GCEM symmetry. 



We point out that there is also a multi-dimensional version of the fluctuation theorem 
where the joint probability distribution of a set of currents (which when summed give the entropy) 
displays a symmetric large deviation function. It is important to note that our symmetry is also different 
from this multi-dimensional case. This becomes clear if the fluctuation theorem obtained in [5] using 
the cycle decomposition approach [12] is considered. Following [8], the fact that the product of the 
transition rates associated with a fundamental cycle in a given direction divided by the product of the 
transition rates in the opposite direction is independent of the states within the cycle, is analogous to 
our condition (|3T|) holding for all columns. 
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Fig. 3 Configuration network of a molecular motor with two possible transition cycles A and B (see text). 



Transitions network with odd Q 

So far we have considered transition networks shown in Fig.[T]with an even number of columns Q. The 
above results can be easily generalized to the case where Q is odd. In this case there are at least two 
adjacent columns with only one configuration. 

In order to do this generalization we can consider an index k that runs only over columns with 
more than one state. Note that wc might also have odd columns with one state. Moreover we take Q 
as the number of columns with more than one state, such that k = 1, . . . , Q. For the case where we 
have an even number of columns and all odd columns have more than one state we have Q — Q/2 and 
k = 2k-\. 

Note that the number of possible complete forward (and backward) cycles is given by Iltf—i n k an< ^ 
only the columns with more than one state are relevant in determining a complete cycle. Now if we 
consider the multiindex i = {ig} and change definition (|32p to 



the formulas (|30|) , (f3"Tj) and (f3"3")) take the same form if Q is odd. Therefore, the proof also works for 
the Q odd case. 



4 Examples 

In order to illustrate how this new symmetry can be established we consider some examples of physical 
systems described by jump processes with a network of states of the type shown in Fig. [TJ where 
the current J r has a clear physical meaning. The examples include a molecular motor in a particular 
network of states, the so-called restricted solid-on-solid (RSOS) model with four sites, and a growth 
model where the nuclcation of the first particle and the completion of the layer take place in the 
time-scales of the same order. 

Molecular motor with four different states 

A molecular motor is a biological protein that converts chemical energy into mechanical work by 
hydrolysis of adenosinctriphosphate (ATP) to adenosinediphosphatc (ADP) and phosphate (P) [53] 
(see also [33] for a review in fluctuation relation for molecular motors). Unlike macroscopic motors, 
which move unidircctionally in a well-defined cycle, a molecular motor performs a random walk driven 
by the chemical potential difference Afi = fi atp — Madp — fJ-p, moving preferentially forward if ATP 
is in excess. This allows one to model molecular motors by stochastic jump processes |26j . Assuming 
that the motor is thermally equilibrated with its local environment the jump rate is proportional to 
exp(/3(Z\/i — W)), where W = Fl is the mechanical work to transport the cargo with force F over the 
distance /. 

The simplest model of a molecular motor would require a cycle of two states. In such a model 
the consumption of ATP would be proportional to the mechanical work. However, realistic molecular 
motors are often characterized by several possible transition cycles, i.e. the motor protein can advance 
by one step through different paths of intermediate configurations with different energy consumption. 
In what follows we consider a hypothetical network of four configurations with two transition paths 



At = 



{A(cf) fc = i,2...,g ; r k = i 




(43) 
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denoted by A and B, as shown in Fig. [3] One interesting question is whether such nano-machine leads 
to a better efficiency at maximum power than a simple linear chain [27] . Here we are interested in 
providing a simple example where our theory can be applied. Clearly, this network belongs to the class 
of networks studied in the previous section with Q = 3 columns. 

If both paths are characterized by different chemical potential differences A\ia and A/jb, the jump 
rates in positive direction will be proportional to Xa = exp [j3 (Aha — Fl)] and Xb = exp [/? (Ahb — Fl)], 
respectively, while the rates for jump in opposite direction do not depend on the chemical potential. 
Moreover, we included a constant factor Y in path B, accounting for a possibly different attempt 
frequency. Finally we assume that the transition rates between cq and c-i are symmetric and equal to 
1, as shown in Fig. [31 The corresponding time evolution operator in a canonical basis {co, ca, cb, C2} 
reads 



C 



/(1 + Xa + YXb) -1 —Y -1 

-X A (I+Xa) -1 

-YX B Y(1 + X B ) -Y 

V -1 -X A -YX b 2 + Y, 



(44) 



We are now going to show that the propagation of the molecular motor defines a time-integrated 
current with a clear physical meaning that exhibits a symmetry in the large deviation regime that 
differs from the one for the entropy. To this end, we note that there are two possible complete cycles 
one going trough ca and the other trough cb- The product of the rates along these cycles in forward 
and backward direction are given by 

Ta = w C(j ^ CA w CA -> C2 w c . 2 ^ Ca = X\ 
Tb = w Co ^ CB w CB ^ C2 w C2 ^ Ca = Y 2 X 2 B 

T B = w Co ^ C2 w C2 ^ CB w CB ^ Ca = Y 2 (45) 



Following the procedure described in the previous section we consider the time-integrated current (|T9l) . 
Note that the polynomials ([25)1 arc of order one for the present model. Comparing the leading order 
terms in Eq. (|29p we obtain the condition 

^MEO) = ^±^. (46) 

Comparing the terms of order a; we obtain a constraint on the transition rates of the form 

Ta + T B Ta^b + T b Xa 



Ta + T b Ta^b+T b \a 



(47) 



where Xa = Xa + 1 and Xb = Y(Xb + 1) are the escape rate in the configurations Ca and Cb, 
respectively. This equation has two solutions, namely Ta/Ta — Tb/Tb, for which ([3"T]) holds so that 
the current is proportional to entropy in the large deviation regime, and = A#, which is the one 
that gives symmetry different from GCEM. Therefore, we need Xa — Xb, which implies 

Y = (l + X A )/(l + X B ). (48) 

Without loss of generality we choose A\ia > A\Ib, which gives Y > 1. In this case, the above restriction 
on Y means that the cycles that go trough configurations with higher chemical potential evolve at a 
slower time-scale. 

Assuming this condition to hold, we consider a time- integrated current of the form (1191) . namely 
the mechanical work 

Jm Fl(^J CA —^ C2 -\- <J CQ ^ CA -(- Jcg->C2 Jco—^cb) (^^) 

for which the factor (|30D is given by 



2FI T A + T B 2FI 



(X 2 A + X 2 )(l + X B ) 2 
(l + X A ) 2 + (1 + X B ) 2 



(50) 
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-0.5 0.5 1 1.5 

Z 

Fig. 4 Scaled cumulant generating functions for F = ft = I = 1, A[la = 3, and A[lb = 2. The black line 
corresponds to the entropy and the red line cor responds to mechanical work (with z —> z/E). They are obtained 
from the minimum eigenvalue of (|51[) and (|52[) . respectively. They are both symmetric and not proportional. 



In order to demonstrate that this quantity exhibits a non-GCEM symmetry, we compute the smallest 
eigenvalues I s (z) and I m (z) of the modified time evolution operators (|T2"j) for both the entropy 



/(l + X A + YX B ) -XX -YX\ 



C s (z) 



-X 
-YX 



(l-z) 



(1-2) 

B 
1 



(1 + x A ) 



-xr z) 





Y(l + X, 



YX 



(l-z) 
B 



-1 

-X'a 
-YXf, 



Y J 



(51) 



and for the mechanical work 



( (l + X A + YX B ) 
-X A cxp(-FZz) 
YX B cxp(-Flz) 



-lexp{Flz) 
(l+X A ) 




V 



-Yexp(Flz) 


Y(1 + X B ) 



-1 



-X A exp(-Flz) —YX B exp(-Flz) 



-lcxp(Flz) 
-Yexp(Flz) 
2 + Y 



(52) 



As expected, the large deviation function for the entropy obeys the GCEM symmetry I s (z) = I s (l — z) 
while the mechanical work fulfills the symmetry I m (z) = I m [E — z) with E given in (|50| . Moreover, 
plotting I s (z) with I m (Ez) in Fig. 0] we see that the rescaled Legendre transforms of the large deviation 
functions arc different. Therefore, the mechanical work is a time-integrated current with a clear physical 
meaning that exhibits a new type of symmetry that differs from the GCEM symmetry of the entropy. 



Restricted solid on solid growth model in a four sites lattice 

The second example is a restricted solid-on-solid (RSOS) model for interface growth 2 Oj - In this 
model the interface configuration is described by height variables hi G Z residing on the sites j of a 
one-dimensional lattice. Particles arc deposited everywhere with rate q while they evaporate with rate 1 
from the edges and with rate p from the interior of plateaus, provided that the restriction hi — hi±i is 
not violated (see left panel of Fig. [S]). On an infinite lattice one observes that the interface roughens 
according to the predictions of the Kardar-Parisi-Zhang universality class [28] . 

Let us now consider a small system with L = 4 sites with periodic boundary conditions. Moreover, 
let us identify all configurations which differ only by translation either in space or in height direction 
so that the interface can be in six possible configurations. By doing that, it can be demonstrated that 
the network of states takes the from shown in the right panel Fig. (see [5D] for details). This network 
belongs to the class of networks studied in the previous section. Moreover, one can easily show that 
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Fig. 5 Restricted solid on solid model. Left: Dynamic rules for (a) deposition and (b) evaporation. Right: 
Transition network of the model with four sites and periodic boundary conditions (see text). 



the escape rates A(c^) = A(c^) = 2q + 2 in column 3 are constant and that the rates across column 1 
obey the condition (|3"7l . 



w 



<■()- 



>C 2 



w 



c - 



>C ^C2— >Cf "'Cf ->-Co 

Hence, as explained in the previous section, the process fulfills the constraint (|33[) . 



(53) 



The most important time-integrated current of the form (|19p with a clear physical meaning is the 
total interface height which increases (decreases) by one whenever a deposition (evaporation) happens. 
The corresponding modified generator is given by 



C h (z) 



{2q + p : 

—qe 
-qe~ 


-2e z 
\ —pe z 



2 -Ape z 
Ap + Aq 


-Aqe~ z 














-e z 
— pe z 
l+p + 3q 

-2qe~ z 
—qer z 



-2qe- z 



~2e z 
-2 + 2q 




-2qe~ 



-2e z 




(54) 



2q ) 



with the factor E 



Mn^ 

4 2p+p 2 



given by relation (|3T)|) . Again its lowest eigenvalue /^(z) has to be 

compared with the lowest eigenvalue I s {z) of the corresponding modified generator for the entropy 
current, which reads 



JC-s(z) 



( 2q + p + 2 


-q z (Ap) x - z 


-q z 





-2q^ z 


-{2qf- z p z 


\ 


-q l - z {Ap) z 


Ap + Aq 





-{Aq) z 
—q z p 1 z 










_ q l-z 





9+1 













-(A q y- z 


— q 1 z p z 


l+p + 3q 


-2q z 


-2 1 - z q z 




-2q z 








-2q l - z 


2 + 2q 







\-(2q) z p 1 - z 








-2 z q x - z 





2 + 2 9 


/ 



(55) 



Plotting the scaled cumulant generating function in Fig. [5] we can see that both are symmetric and not 
proportional. Unfortunately, for L > 4 sites the height current is no longer symmetric because in this 
case the network of transitions is not of the form shown in Fig. [T] In the next example, we introduce 
a model where the space of states is of the form displayed in Fig. [T]for any system size. 



Growth process with instantaneous monolayer completion 

As a third example, let us define a growth process with deposition and evaporation which has the 
special property that after a nucleation the actual monolayer is completed on a very short time scale. 
The model is defined on a one-dimensional lattice with particles of species a and /3, assuming that 
only a particles can be deposited on top of /3 particles and vice versa (see Fig. [7]). Once a particle is 
deposited with rate e?" at site i on a flat surface, the subsequent event is either the completion of the 
layer or the evaporation of the particle with rate ef'^. This describes a limit in which a monolayer is 
completed almost instantaneously after the first deposition. The transition rates for the completion of 
a layer is 1 and the transition rate for the reversed event, that is, the evaporation of L — 1 particles of 
a monolayer, is e. 
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200 




Fig. 6 Scaled cumulant generating functions for the RSOS model for p = 0.02 and q — 10. The red line 
corresponds t o he ight ( with z — > z/E) and the black line to the entropy, they are obtained from the minimum 
eigenvalue of (|54|l and (|55|l . respectively. 



_□ 









3 


a 


a 


a 




Z 













Fig. 7 Possible cycle of transitions for the case L = 4. The a particles are blue and the j3 particles red. The 
left-most configuration is a flat interface with /3 particles. 



The network of configurations of this model is of the form shown in Fig. [T] with Q = 4, where 
Co and C2 correspond to flat interface configurations while the columns c\ and C3 have L states, each 
corresponding to a single particle nucleated at one of the L sites. In this model the physically relevant 
time-integrated current of the form (|T9)) is the interface height, which increases by 1/L when a particle 
is deposited on a flat interface, and by (L — 1)/L whenever a monolayer is completed. 



The product of rates on a forward and backward cycle through c\ and are 



df£ and 

1 J 



while the escape rates are X(c\) = Xf = 1 



and A(c|) = A£ = 1 + e£. Hence, the 




Fig. 8 Scaled cumulant generating functions for the growth model with fast monolayer completion. We used 



L = 3, df = 1, df = 2, df = 3, d? = 2, d£ = 2, dg = 5, ef = ef = ef 
The black line is related to entropy and the red line to height. 



3, and e = 0.05. 
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polynomials (|26| with degree 2L — 2 are given by 

L L 



/(*) = EE^IK a 



i) 



i) 



i=l 3 = 1 



and for the factor (|30[) we get 



In 



E 



d a d p 



1,3=1 e i e i 



The constraints (f3"3"| take the form 



E 



i,j=l a i a j 



i,j=l a i a j G k A\ , . . • , A^_i: A i+1 , . . . , A L , Aj . 



' A J-1' A J + 1' 



,2 V L P ce p' 3 2 a /3 



A, 



A i-1' A i+1' • • • ) A £j A x , . . . , Aj_ t , A J + 1 , 



* ' L 



' L 



(56) 
(57) 

(58) 



(59) 



where k = 1, 2, . . . , 2L — 2. Whenever the deposition and evaporation rates fulfill the above equation, 
the probability distribution of velocity will symmetric with respect to (|58[) . Again, this symmetry is 
generally different from the GCEM symmetry. As an example in Fig. [5] we plot the scaled cumulant 
generating function for the entropy and height for a system with L = 3 sites, where the transition rates 
are chosen in such a way that the escape rate is constant in columns c\ and C3 so that the relation 
([BUI) is satisfied. 



5 On the origin of the symmetry 

In the following we demonstrate that the new symmetry we proved here also comes from time-reversal. 
However, in order to see the symmetry we have to consider a functional of a group of trajectories 
instead of a single trajectory. We will restrict our analysis to the four-state system shown in Fig. [31 but 
we expect the same kind of proof to be valid for general networks of the form shown in Fig . [TJ Before 
going into detail, we provide a simple demonstration of the GCEM symmetry, where it becomes clear 
that its physical origin is related to time-reversal. 

A simple demonstration of the fluctuation theorem 

As in Sect. 2 we consider stochastic trajectory with M jumps in the time interval [to,tf] and jumps 
from c(ti) to c(tj-|-i) take place at times £$. We denote this trajectory by M t and its weight is given 

by 

M-l 

W0 M ,t] = exp[-X(c(t M ))(t f - t M )\ J] w c{u) ^ c{u+l) exp[-X(c(ti))(t i+1 - t,)], (60) 

i=0 

where A(c) is the scape rate from state c. We note that we should also multiply the weight by the 
probability distribution of the initial state, but we are considering a uniform initial distribution of 
states. The reversed trajectory, where the system starts at state c(^m) at time i/ (again with the 
uniform distribution) and jumps from state c(U+i) to state c(U) at time tf + to — U+i, is denoted by 
It M,t f +t -t- The weight of this trajectory is written as 

A/-1 

W[tJ M,t f +t -t] = exp[-X(c(t M ))(tf - t M )] \\ w £ (i l+1 )-> c (i,) exp[-A(c(^))(t i+ i - U)]. (61) 

i=0 
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The entropy current (J7]) is related to the weight of a trajectory divided by the weight of the time- 
reversed trajectory by 

i t r7^ m W\yM,tf+t -t] ,„_,. 

exp(-J s [C M ,t}}) = — T : , (62) 

W[C M ,t\ 

From dB2J it follows that 

W0 M ,t] exp(-X)5(J a [d Mtt ] -X) = W[t M,t f +t o -t]S{J s M,t] - X) . (63) 

Summing over all possible trajectories and using J s [ft M,t] = — Js\& M,tf+t -t] we obtain the fluctuation 
theorem ((5J), i.e., 

%=^ = ^l-X). (64) 

This demonstrates that the GCEM symmetry is a direct consequence of the fact that the entropy is 
the weight of a trajectory divided by the weight of the time-reversed trajectory. 



Time-reversal of a group of trajectories 

We now consider Markov jump process with the network of states given in Fig. [31 In the following we 
will denote the states ca and Cb by cf and cf , respectively. Here, instead of considering one stochastic 
trajectory, we consider a certain group of trajectories (or class of trajectories), which can be defined 
as follows. Two trajectories belong to the same class if they have the same number of jumps taking 
place at the same times ti, following of the same sequence of columns. For example, the trajectories 
Co — > cf — > C2 — > Co and Co — >• cf — s- C2 — > Co belong to the same class. More generally, a trajectory 
that goes trough column 1 for K different times pertains to a class with 2 K trajectories. 

We denote this group of trajectories by {C M,t] an d its weight, which is the sum of the weights of 
each trajectory in the group, by R[{ Cmj}]- We define a quantity J such that 

i t\(V$ n\ M,t f +t -t}} , > 

exp(-J[{C M ,t}]) = mr ^ J , (65) 

n[{U M,tj\ 

where R[{U M,t f +t -t}] is the sum of the weights of the reversed trajectories in the group and are using 
a tilde to denote functionals of the group of trajectories. Note that, unlike the entropy, the current J 
is a functional of a group of trajectories. What we show next is that in the case Ai = A(c^) = A(cf ) 
(which is the sufficient condition (|34[) for the symmetry), then the the functional J r , which is generated 
by the current J r in the way explained below, is equal to J, provided the increments Ok defined in (|19|1 
are chosen appropriately. 

The difference between a current that is a functional of the stochastic path ~C M,t and a current 
that is a functional of the group of paths {C M,t} is i n the increment when there is a jump to an state 
at column 1. For the second kind of currents the increment is as follows. Let us consider a trajectory 
which at time tj_i hops to a state c(Aj) that pertains to the column ci, and stays in that state during 
the time interval ti — t%—\. The contribution to the current </[{C^A/.t}] is given by 

™c(t i - 1 )^c? w c?-+c(t i+1 ) exp[-A(c^)(tj - + w c(U i) ^ c b w c b _). c ( ti+1 ) exp[-A(cf )(U - t,_i)] 
w c(t i+1 )^cf ? V- ! -c(t i _ 1 ) exp[-A(cf ){U - ti-i)] + w c{t . +l) ^ c BW c B^ c{t ._ i) exp[-A(cf )(U - U-i)] ' 

(66) 

If the escape rate is constant at column 1 the above term becomes 

w c{t i+l )^ct W ct^c(U^) + w c(U +1 )^cf W c? _, c(ti _ l} ' 

Note that, the same relation holds if the relation ([3T)) is valid at column 1. The fact that the scape 
rates are constant simplify the situation considerably: in this case the dependence on the time interval 
ti — tj_i disappears. From now on we assume that the escape rates are constant at column 1. 
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The current J r is a functional that is invariant within the class of trajectories, i.e., all the trajectories 
in the same group have the same value of J r . This important property comes from the fact the the 
current (|19[) is defined in a way such that it does not discriminate between different states in the same 
column. Hence, J r induces a current J r defined on a class of trajectories by the relation 

Jr\{C M,t}] = Jr[^M,t) for ~d M ,t G M ,t}- (68) 
If we consider the current J r with the increments 

00 = 0! = -In c °^ Cl Cl ^° 2 C0 ^ Cl Cl ^ C2 , 2 =\ n ^^, (69) 

then it is clear that J r [{C^M,t}] is equal to the current given in (p5|) . It then follows that 

R[0 M ,t}} exp(-X)5(j r [0M,t}} ~ x ) = R[0M,t f +t o -t}Wr[0M,t}} - x ) • (70) 

Now summing over all possible group of trajectories we get P if — cxp(— X). Since relation (|68p 

gives P(J r = X) = P( J r = X) we finally obtain 

P( J r = -X) , 

This last relation implies the symmetry for the current J r with the choice of increments ((69"]) . However, 
as we proved in (|22j) . at large time all the currents J r are proportional for any choice of increments 0k- 
Therefore, this argument shows that the symmetry in the current J r comes also from time-reversal, but 
at a more coarse grained level. We point out that functionals of a group of paths have been considered 
also in [29], however in a rather different context. 

We believe that the same demonstration can be extended for the more general network of Fig. [T] 
However, formalizing the proof for this more general case leads to much more cumbersome formulas. 
We pretend to address this problem in future work. 



6 Conclusion 



In this paper we have investigated a new symmetry of time-integrated currents which is generally 
different from the GCEM-symmetry of the entropy current in the long time limit. It is valid for a 
very restricted class of Markov processes in the sense that they have the very peculiar network of 
configurations shown in Fig. [T] Moreover, the symmetry appears only when the transition rates fulfill 
a certain set of constraints. Nevertheless, to our knowledge, this is the only case where a current not 
proportional to the entropy in the large deviation regime displays a symmetric large deviation function. 

We showed three physical examples where this current is a relevant physical observable. We consid- 
ered a toy model for a molecular model, where the mechanical work has a symmetric large deviation 
function different from the one for the entropy. Moreover, we analyzed two growth models where the 
height displays such a symmetry. 

As is the case of the GCEM symmetry, the origin of our symmetry seems to be related to time- 
reversal. However, as we showed in Sec 5 for a 4 states system, the symmetry becomes clear only 
when we consider a functional of a group of stochastic trajectories. That is, the symmetric current 
as a functional of the group of trajectories is given by the weight of the trajectories divided by the 
weight of the reversed trajectories. We demonstrated the symmetry with this grouping of trajectories 
argument only for the specific 4 states system, but we expect the same kind of proof also to be valid 
for the whole class of networks considered here. 

A natural extension of the this work is to look for other currents, in a more general space of states, 
with a symmetric large deviation function. This can be done by looking for the conditions that the 
transition rates have to fulfill in order for the characteristic polynomial of the modified generator 
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to be symmetric. We expect that whenever this characteristic polynomial is symmetric, the origin 
of the symmetry should be related to time-reversal (of a trajectory or a group of trajectories). One 
might speculate about a situation where the characteristic polynomial of the modified generator is not 
symmetric, however its minimum eigenvalue is still symmetric. In this case the origin of the symmetry 
might be also related to time-reversal, but of some most probable stochastic trajectory dominating a 
sum over different trajectories. Ultimately, it would be of great theoretical interest for nonequilibrium 
statistical physics to find out which are the time-integrated currents with a symmetric large deviation 
function and what might be the origin of these symmetries. 
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